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ABSTRACT 

We present empirical relations describing excess emission from evolved stars in the 
Large Magellanic Cloud (LMC) using data from the Spitzer Space Telescope SAGE 
(Surveying the Agents of a Galaxy's Evolution) survey which includes the IRAC 3.6, 
4.5, 5.8 and 8.0 /um and MIPS 24, 70 and 160 /um bands. We combine the SAGE data 
with the Two Micron All Sky Survey (2MASS; J, H and Kg) and the optical Magellanic 
Cloud Photometric Survey (MCPS; U, B, V and I) point source catalogs in order to 
create complete spectral energy distributions (SEDs) of the asymptotic giant branch 
(AGB) star candidates in the LMC. AGB star outflows are among the main producers 
of dust in a galaxy, and this mass loss results in an excess in the fluxes observed in 
the 8 and 24 fj,m bands. The aim of this work is to investigate the mass- loss return by 
AGB stars to the interstellar medium of the LMC by studying the dependence of the 
infrared excess flux on the total luminosity. We identify oxygen-rich, carbon-rich and 
extreme AGB star populations in our sample based on their 2MASS and IRAC colors. 
The SEDs of oxygen- and carbon-rich AGB stars are compared with appropriate stellar 
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photosphere models to obtain the excess flux in all the IRAC bands and the MIPS 24 
Hm band. Extreme AGB stars are dominated by circumstellar emission at 8 and 24 
lim thus we approximate their excesses with the flux observed in these bands. We 
find about 16,000 0-rich, 6300 C-rich and 1000 extreme sources with reliable 8 //m 
excesses, and about 4500 0-rich, 5300 C-rich and 960 extreme sources with reliable 24 
fim excesses. The excesses are in the range 0.1 mJy - 5 Jy. The 8 and 24 /xm excesses for 
all three types of AGB candidates show a general increasing trend with luminosity. The 
color temperature of the circumstellar dust derived from the ratio of the 8 and 24 fim 
excesses decreases with an increase in excess, while the 24 fim optical depth increases 
with excess. The extreme AGB candidates are the major contributors to the mass loss, 
and we estimate the total AGB mass-loss return to the LMC to be (5.9 — 13) x 10"'^ 
Mq yr-i. 

Subject headings: (galaxies:) Magellanic Clouds, infrared: stars, stars: AGB, stars: 
carbon, stars: mass loss 



Introduction 



Towards the final stages of their evolution, low- and intermediate-mass stars (0.8 - 8 Mq) enter the 
asymptotic giant branch (AGB) phase. Characterized by cool (~3000 K) photospheric temperatures 
and dusty circumstellar (CS) shells which absorb optical light from the star and re-emit it at longer 
wavelengths, AGB stars are among the brightest infrared (IR) sources in a galaxy. During the 
AGB phase, a significant fra ction of the star's original mass is lost at rates of up to -^10"^ Mq yr~^ 
(|Wood. Bessell &: Foxl Il983l : ?), resulting in the chemical enrichment of the interstellar medium 
(ISM). 

AGB mass loss is believed to be a two-step process: pulsations first levitate material above the 
photosphere, where the cool temperatures result in the formation of dust grains. Radiation pressure 
then drives t he dust grains (which are collisionally c oupled with the gas) outward in an efficient 



stellar wind (IWickramasinghe. Donn &: Stecheiill966l : iGoldreich Sz Scovilldll976l : iBowen Sz Willson 
199ll : IWachter et al.lbooi T An increase in luminosity (and hence radiation pressure) must therefore. 
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The formation of a circumstellar envelope (CSE) is a direct result of the mass-loss process. 
The flux from the CS dust shell appears as mid-infrared (MIR) emission in excess of that due 
to the central star alone. This MIR excess is directly related to the rate of mass loss, and is 
therefore expected to increase with increasing luminosity of the central star. Various studies in 
the past have demonstrated the relationship between the MIR excess and the rate of mass loss. 
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Skinner &: Whitmord (|l988l ) showed that the MLR derived from CO rotational transitions was pro 



Jura 



portio nal to the strength of the 9.5 //m siUcate dust feature in opticahy thin 0-rich stars 
(|l987l ) derived a rel ation between the MLR and 60 //in exces s for O- and C-rich AGB stars in the 
solar neighborhood. iKnapp. Gunn fc Wvnn- Williams! (119921') found tha t the 1 2 //m emission from 
elliptical galaxies was proportional to the measured MLR. I At he v et al.l (|2002l ) compared the MIR 
excess emission from 9 galaxies to that of Galactic and LMC AGB stars and derived a proportion- 
ality relation between the MIR excess and MLR. 

While the study of Galactic AGB stars is hindered by obscuring Galactic dust along the line 
of sight and high uncertainty in distance estimates, the Large M agellanic Cloud (LMC) is ideal for 
such observations owing to its proximity (~50 kp c. lFeastlll999l) and favorable viewing angle (35°, 



van der Marel k Cioni"2001 



1981 



. Early LMC survevs (IWe sterlund et al 



1978: 



Westerlund. Olander &: Hedin 



Eebeirot et al...l983l : iBlanco. Blanco fc McCarthy. 1980 : Blanco fc McCarthvlil983, : .Frogel fc Blanco 
19901 ) looked for AGB stars at optical and near-infrared (NIR) wavelengths, and hence preferen- 
tially detected sources w ith optically thin C S dust shells. The Infrare d Astronomy Satellite (IRAS; 
Neugebauer et al.lll984l ) survey of the LMC (ISchwering &: Israellll990l ) in conjunction with ground- 



losing evolved stars (jLoup et al 



1997 



based NIR confirma tion (IReid. Tinnev &: Mouldlll990l) he l ped identify several of the brightest mass- 



Ziilstra et al.lll996l ). iTrams et al.l (jl999l ) performed follow-up 



mid-infrar ed (MIR) photomet ry and spectroscopy of 57 sources using the Infrared Space Observa- 
tory (ISO: iKessler et al.lll996l ) in order to determine the chemical composition of the CSEs, and 



van Loon et al 



(|l999l ) used radiative transfer modeling to derive luminosities and MLRs for these 
sources. More recently, the LMC has been s urveyed in the optical by the Magellanic Clouds P hoto- 
metric S urvey (MCPS: Izaritskv et al.l ll997l'). and in the NIR by DENIS (lEpchtein et al.lll994ll and 
2MAS S dskrutskie et ahlboodl. T he LM C data from these two surveys can be found in lcioni et all 



performed bv lEgan. van Dvk &: Pried (|200ll ). 



(|2006l ) and iNikolaev &: Weinberg! (120001). re spectivelv. A mid-infrared survey using the Midcourse 
Space Experi ment (MSX: !Price et al.!!200l!) . which was four times more sensitive than IRAS, was 



As part of the SAGE survey (Surveying the Agents of a Galaxy's Evolution; !Meixner et al 



2006 ). a ~7° X 7° area of the LMC was imaged in the Spitzer Space Telescope ( Spitzer AWernei et al 
20M) IRAC (3.6, 4.5, 5.8 and 8.0 fxm) and MIPS (24, 70 and 160 //m) bands. One of the main goals 
of SAGE was to detect all the evolved stars with MLRs > 10"^ Mq yr~^, to characterize the total 
rate at which material is returned to the ISM by dusty evolved stars, and to better understand the 
physics governing mass loss among evolved stars in the LMC. When complete, the SAGE data will 
be ~1000 times more sensitive than MSX, and they will allow a detailed quantitative derivation of 



MIR spectroscopic observations of evolve d stars (see, e 
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2006! ). Continuing the analysis begun bv !Meixner et al.! (|2006! ) . !Blum et al.! (|2006! ) identified about 



van Loon et al. 


1999, 
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32,000 color-selected evolved stars brighter than the tip of the red giant branch (11.85 mag in the 
IRAC 3.6 fim band), including 17,500 oxygen-rich (0-rich), 7000 carbon-rich (C-rich), and 1200 
"extreme" AGB stars, and presented color-magnitude diagrams (CMDs) of the SAGE epoch 1 data. 
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Our aim is to investigate the MLRs of the AGB candidates selected from the SAGE survey, 
which will require detailed radiative transfer modeling of the CSE around each such star. In order 
to simplify this effort to model ~ 10^ stars, we would like to constrain the range of some of the 
input model parameters, such as the total source luminosity, the typical temperature of the CS 
dust and the optical depth at a given wavelength. As a first step towards this modeling goal, 
therefore, in this paper we take an empirical approach and calculate MIR excess fluxes due to dust 
from these stars. The distance independer ice of the MIR colors o f AGB sources makes them a good 
choice for tracing the mass loss (see, e.g.. IWhitelock et al.ll2006l ). but our study based on the MIR 
excesses for sources in the LMC is aided by the fact that the distance to the stars in the LMC is 
essentially the same. Data from the SAGE survey also provides an unprecedented opportunity to 
include a large number of AGB star candidates. We will show that the MIR excess can be used as 
a proxy for the MLR. In a very statistical way, we study the overall trends of excess and derived 
quantities such as color temperature of the dust and MIR optical depth with source luminosity. 
In subsequent papers, we will also be presenting the luminosity function and detailed radiative 
transfer calculations for the AGB star candidates selected in this paper. The paper is organized 
as follows: in ^we describe the SAGE database and our observational sample of AGB stars. Our 
procedure for calculating IR excesses is explained in ^j3l We present our results in ^ and in ^ we 
compare our results to previous work and discuss their implications for future AGB studies in the 
LMC. 



2. Data 



The SAGE catalog and archive point source lists from both epochs of observations have been deliv- 
ered to the Spitzer Science Center (SSC), a nd are available for do wnlo ad at the SScEI. Th e SAGE 
epoch 1 point source catalog is discussed by lMeixner et al.l (j2006l ) and iBlum et al.l (j2006l ). In this 
study, we select sources from the IRAC epoch 1 archive instead of the catalog. The archive accepts 
fainter sources than the catalog. Faint limits for both epochs are 18.5, 17.5, 15, and 14.5 mag for 
IRAC 3.6, 4.5, 5.8, and 8.0 /xm respectively, compared to 18, 17, 15, and 14 for the catalog. A 
source is excluded {i.e., culled) from the archive if it has neighbors within a 0.5" radius, whereas 
this radius is 2" for the catalog. This cullin g procedure erisures the creation of reliable lists of point 
sources at the expense of completeness (see lMeixner et al.ll2006l ). In order for two sources to appear 
blended, they both have to be detected, and they have to be bright enough. The probability of such 
blending in the 24 /xm band is very low, since the 24 /im catalog has very few sources (~ 40, 000) 
compared to the IRAC catalog. However, we do find a few (~40 out of ^ 10^) IRAC sources that 
are matched to more than one MIPS 24 fim source which we corrected in our results and in the 
main catalog. In addition to this culling, a flux value for each archive source in any of the IRAC 
bands is non-null only if its signal-to-noise ratio (S/N) is greater than 5. For the catalog, the S/N 



http://ssc.spitzer.caltech.edu/legacy/sagehistory.html 
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is 6 for the IRAC [3.6], [4.5], and [5.8] fim bands, and 10 for the [8.0] fim band. As a result of 
these criteria, the archive has shghtly more sources and more flux values than the catalog, allowing 
the inclusion of more faint AGB candidates. As of the second delivery, version S13, the archive 
has ~4.5 million sources compared to ~4.3 million in the catalog. See the SAGE data delivery 
documenll^ for details of the source selection in the catalog and archive. The fact that the fluxes 
of these fainter AGB candidates are more uncertain when using the archive list is mitigated by our 
requirement of 2MASS detections thus affirming that the point source is real. The AGB source 
lists are extracted from a universal table which lists the nearest neighbor in the MCPS and IRAC 
epoch 1 archive to each source in the MIPS24 epoch 1 full catalog ('^40,000 sources). The matching 
radius is 3", but the nearest neighbor is used in the match. The IRAC archive sources are also 
bandmerged with the 2MASS catalog. 



We classify sources as low- or moderatel y-obscured 0-rich/C-rich AGB can didates based on 
their location in the K„ versus (jCioni et al.l l2006l : I Blum et al.l l2006l ) , but we exclude 

stars without H-band fluxes from our list, as our procedure for calculating the excesses (see ^ relies 
on an H-band detection. This criterion probably excludes the faintest 0-rich AGB stars, and thus 
does not significantly affect our results. As a separation into 0~rich a nd C-rich chemistries based o n 
near-IR colors is not possible for the heavily obscured fJ-[3.6] >3.1. iNikolaev &: Weinbergl ([2000)) 
"extreme" AGB candidates, we follow a procedure similar to iBlum et al.l (j2006l ) and select these 
sources based on their location in the [3.6] versus J-[3.6] CMD or (when no 2MASS counterpart 
exists) the [8.0] versus [3.6]-[8.0] CMD. A simple trapezoidal integration of the optical U through 
MIPS24 fluxes is performed to estimate the luminosities for all our sources. 



The lKastner et al.l (|2008l ) list of objects wi th spectroscopi c class ifications contains 14 sources 
classified in our study as 0-rich AGB stars. iKastner et al.l (j2008l ) classify 10 of these as red 
supergiants (RSGs). We also find that we have classified 21 of their sources as extreme AGB stars 
- most of these have been identified as possible HII regions, none of them are RSGs. Since these 
are point sources, we interpret HII regions as compact HII regions or massive young stellar objects 
(YSOs). We thus realize that our sample of AGB stars may be contaminated by non-AGB objects 
such as RSGs and YSOs. We estimate this contamination using simple color-magnitude cuts. 



The figures in lBlum et al.l (|2006l ) also show that the AGB stars and RSGs are not well separated 
in the MIR CMDs. It is possible that some of the most luminous sources in our list are RSGs. 
We find 556 objects (124 classified as 0~rich, 17 as C~rich, and 415 as extreme AGB candidates) 
more luminous than Mboi = —7.1, the classical AGB luminosity limit. It is not unlikely that some 
of these sources are AGB stars. Luminosities above the classical limit can be achieved by AGB 
stars at the peak of their pulsation cycle (as in the case of OH/IR stars) or by stars undergoing 



^http://data.spitzer. caltech.edu/popular/sage/20080204_enhanced/documents/SAGEDataDescription_Delivery2.pdf 
^The UBVI and JHKs magnitudes are dereddened using the same procedure as in Cioni et al. ( 20061 ) 
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hot bottom burning (HBB. iBoothrovd &: Sackmannlll992l ). Nevertheless, these numbers provide a 
very conservative upper hmit for the RSG contamination in our sample. 

The redde s t sour ces in our list fall in the region of the MIR CMD also populated by YSOs. 
Whitney et alJ (120081 ) isolate regions in the [8.0] vs [8.0]-[24] CMD occupied more densely by 
YSO models (Figure 3 in their paper). This separation includes a stringent cut at [8.0]-[24]~2.2 
(corresponding approximately to a 24 fim flux > the 8 /im flux) to exclud e AGE stars. We find that 
5 0-rich and 24 extreme AGB star candidates in our list are part of the IWhitnev et al.l (|2008l ) list 
of high-probability YSOs. A more conservative estimate is obtained by looking for sources fainter 
than [8.0]=7 with F24f_im / Fsum ^ 1, which puts the YSO contamination in our lists at about a 100 
sources (61 0-rich, 14 C-rich, and 34 extreme AGB candidates). 

Figure [1] shows the Kg versus J-Kg Hess diagram for the region populat ed by AGB stars. The 
IRAC and MIPS24 CMDs are shown in Figures [2] and El IBlum et al.l pOod ) noted the presence of 
a fainter, redder population of 0-rich sources (finger "F" in their Figure 6). The same population 
can be seen in our IRAC-MIPS24 CMD (Figure [3]) at [8.0] magnitudes fainter than ~10. We 
use a magnitude cut at 10.2, shown in the figure as a solid line. Throughout this paper, we will 
differentiate between the bright and faint 0-rich populations based on this magnitude cut. Almost 
80% of the 0-rich stars in our sample are fainter than [8.0] =10.2. Photometric information for a few 
sources of each type is shown in Table [H The electronic version of this table contains information 
for our entire list of AGB star candidates. 

Figure H] shows the luminosity function for all three classes of AGB candidates. The lower 
luminosity limit {Mf,oi ~ —4 ) for the 0-rich and C-rich AGB candidates is set by our color- 
magnitude cuts which exclude sources fainter than the tip of the RGB. The C-rich sources have 
a tighter luminosity distribution than their 0-rich counterparts (there are only a handful of C- 
rich sources brighter than Mboi = —7.5, whereas the 0-rich distribution falls off around Mhd = 
—8.5). The range in luminosities for the extreme AGB candidates is M^oi = —4.5 to Mboi ~ —10. 
The extreme AGB candidates thus have the highest luminosities in the sample. However, we 
have insufficient information at this point to say anything concrete about the breakdown of these 
sources into 0-rich and C-rich chemistries. The spectroscopic follow-up to SAGE will provide some 
information about the dust chemistry of these extreme AGB candidates. These distributions peak 
at Mhoi « —5 (0-rich), M^oi ~ —6 (C-rich) and M^oi « —6.5 (Extreme). The vertical dashed line 
shows the c lassical AGB luin inosity limit. While this limit can be exceeded by deeply embedded 
AGB stars (jWood et al.lll992l ). 0-rich sources with luminosities higher than Mhoi = —7.8 are more 
likely to be RSGs, while the brightest extreme AGB stars may be massive YSOs. We will discuss 
the astrophysical implications of these luminosity functions in a future paper. 
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3. Procedure 

We estimate the MIR excess emission of the low- and moderately-obscured 0-rich and C-rich AGB 
star candidate populations by comparing their observed SEDs to an expected SED for the stellar 
photosphere. The IR excess in each wavelength band is calculated by comparing the total flux from 
the source observed in that band to the flux expected from the central star as prescribed by a "best- 
fit" model photosphere. As the emission from CS dust dominates the MIR flux of extreme AGB 
stars, we set the MIR flux equal to the excess in each band for our extreme AGB star candidates. 
In this work, we are interested in describing overall trends in the AGB parameters such as IR excess 
as opposed to detailed models of each source. To this end we choose a single model photosphere 
for each type of AGB star - one that best fits the S EP shape of AGB stars wit h little or no dust. 



We use the plane-parallel G -rich MARCS models of iGautschv-Loidl et al.l (j2004l ) and the spherical 



0-rich PHOENIX models of lHauschildt. Allard k BaronI (|l999l ) to calculate the photospheric AGB 



star emission. The differences between plane-parallel and spherical models are insignificant as far 
as the resulting broad band AGB star photospheric information is concerned. The range of model 
parameters we search for a "best-fit" are as follows: the thirty-two solar-metallicity C-rich AGB 
models have C/0 ratios between 1.1 and 1.8, effective temperatures ranging from 2600 K to 3200 
K, and surface gravity logg between about —0.76 and 0. The ~200 0-rich AGB models of solar 
mass and solar metallicity have effective temperatures ranging from 2000 K to 4700 K, and logg 
between —0.5 and 2.5 in steps of 0.5. Synthetic photometry in the optical, near- and mid-infrared 
was obtained from each of these models by convolving their SEDs with the filter response curves of 
the Johnson-Kron-Cousins UBViS , 2MASS JHKsl , and Spitzer IRAC0and MIPS 24 ^m|l| filters. 
The oxygen-rich models had spectral information in the range 0.1-1000 /im and the convolutions 
were done directly with the models. In order to compensate for the insufficient wavelength coverage 
(~0.5-25 fim) offered by the carbon-rich models, the U and B band fiuxes were dropped, and the 
flux in the MIPS24 band was extrapolated to 30 fim assuming a Rayleigh-Jeans falloff. 

As we are most interested in obtaining the correct shape of the SED for a best-fit, the model 
SEDs are scaled to the H-band fiux of the median SED of the hundred bluest sources in V— Kg 
color. We use the same model to calculate the excesses for both the bright and faint 0-rich 



The MCPS magnitudes were placed on the Johnson-Kron-Cousins UBVI system. 



The detector quantum effi- 
ciency curve was obtained from the Las Campanas Observatory website, http://www.lco.cl/lco/index.html Filter 
profiles for the Johnson U, H arris B, V, and Cousins I filters were obtained from the references in Table 9 of 



Fukugita. Shimasaku fc Ichikawa ( 199 5). 



'^The 2MASS fiher relative spectral 
tained from the 2MASS All-Sky Data 



Cohen et al 



responses (RSRs) derived by 
Release Explanatory Supplement on the 



1 20031 ) were ob- 
worldwide web at 



http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec64a.htmli 



The IRAC RSRs are plotted in iFazio et al.' ('20041), and were obtained from the Spitzer Science Center IRAC 
pages at http:/ /ssc. spitzer. caltech.edu/irac/spcctralresponse. html 

^The MIPS dRieke et all 



I2OO4I) 



RSRs were obtained from the Spitzer Science Center MIPS pages at 



http://ssc.spitzer.caltech.edu/mips/spectralresponse.html 
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populations because the bluest 0-rich candidates lack detections in the 24 ^um band and thus can 
not be separated into bright and faint populations. The model that comes closest to describing 
the oxygen-rich median SED has surface gravity log 5=0 and an effective temperature r=4000 K. 
The corresponding carbon-rich best-fit model parameters (for a C/0 ratio of 1.3) are log ^=-0. 43 
and T=3200 K. Figure [5] demonstrates the filter-folding for these best-fit models (top), and shows 
the comparison to the median SED of the bluest sources (bottom). There are many more 0-rich 
candidates that are fainter and bluer than most of the C-rich candidates, which might explain the 
considerable difference in the best-fit model temperature and radius of the two types of sources. 

To determine the excess flux Xj^ in a band centered around frequency v due to CS emission, 
the best-fit model flux must first be scaled to the flux of each source at some "pivot" wavelength 
and then the difference between the observed flux F^^ and the corresponding scaled model flux in 
that band must be calculated. Fitting the SEDs of all our sources to one best-fit model photosphere 
is a simple first approach towards modeling the CSE around each AGB star candidate. The effects 
of interstellar and CS extinction on our sources is minimal in the 2MASS JHKg bands, which also 
contain the wavelength range corresponding to maximum emission from AGB star photospheres !§ 
Almost all of our 0-rich sources peak in the H bancil, while over two-thirds of the C-rich sources 
peak in the Kg band. All of the photospheric models available to us exhibit H-band SED maxima. 
For these reasons, we place the pivot wavelength in the H band (centered at 1.65 /xm). The excess 
then depends on the source flux Fy and model flux F™"*^ according to 



where the subscript H denotes H-band fluxes. This equation may overestimate the MIR excesses 
from the redder C-rich sources with Fk > Fh- On the other hand, using the same best-fit model 
with a Kg-band pivot would produce underestimates for the excesses. For the reddest C-rich source 
in our sample, the difference in 8 //m excess resulting from choice of pivot is ~ 30%. This is to be 
regarded as an upper bound to the error introduced in the excess determination due to the choice 
of the H-band as pivot. While this will not alter the general trends we discuss in this paper, it will 
affect our numerical results for higher excesses (mass-loss rates). 

The photometric errors associated with the source fiuxes were used to estimate the error 6X 
in the calculated excess. An excess measurement in a wavelength band centered around frequency 
V was deemed "reliable" only if its signal-to-noise ratio was greater than 3. In other words. 



Figure [6] illustrates the effect of selecting sources using this "3-a"" criterion. The distribution 
of sources that are rejected based on this criterion is fairly symmetric about zero, with a slight 



The most obscured sources will suffer from CS extinction even in the NIR bands, but choice of pivot wavelength 
is not an issue for these sources, as they are probably members of our extreme AGB list. 

^This is partly due to the fact that the opacity of the H~ ion reaches a minimum in the H-band. 
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asymmetry on the positive excess side, indicating that our cut is conservative. A similar distribution 
of rejected sources is seen for the C-rich candidates, but they are substantially fewer in number. 
In both cases, sources with excesses below ~ 0.1 mJy are rejected. 

Very little can be inferred about the chemical composition of the dust shells around the extreme 
sources without spectroscopic confirmation, although most of these objects are pr obably C-rich 



We find that 60 of our extreme AGB candidates are identified as C~rich in the Kastner et al 



(j2008l ) study, while only 9 are classified as 0-rich. Our procedure for calculating excesses relies on 
a classification into 0-rich or C-rich sources, which is not possible for most of the stars in this list. 
However, the excess emission from their extremely dusty CSEs dominates over the photospheric 
emission in the MIR, so that we can set the excess equal to the MIR fiux to a good level of 
approximation. We find ~ 8200 0-rich, 5800 C-rich and 1400 extreme sources with reliable 8 
/im excesses, and about 4700 0-rich, 4900 C-rich and 1300 extreme sources with reliable 24 fim 
excesses. Table [2] shows the excesses calculated for the fifteen sample sources shown in Table [H 
The electronic version of Table [2] lists all the sources with valid 8 fim excesses. 

We estimate the temperatu re of the CS dust and th e opt ical depth fr o m the 8 fiui and 24 fim 
excesses in a manner similar to iThompson et al.l (j2006l ) and iDaval et al.l (|l998l ) . The continuum 
dust emission is modeled as a blackbody at temperature T^, with optical depth t^, 



where 



TT 



R 



IT 



(3) 



(4) 



is the geometrical dilution factor with the distance from the central star corresponding to 
maximum emission due to dust at both 8 and 24 /xm. For small optical depths, the excess is 
proportional to the opacity k,^. If the emissivity of the dust can be modeled by a power law in 
the relevant range of wavelengths, we have 



(5) 



Assuming a power-law emissivity in the MIR ignores any effects due to strong absorption or 
emission features from silicate dust for example, which can be significant for the more obscured 
0-rich source s. The power-law index will also depend on the dust species in general (see, e.g., 



Hofner 1120071 ). The A ^ power-law dependence adequately describes the emissivity at MIR wave- 



lengths for the purpose of this paper (the Spitzer 8 and 24 /im bands will only detect the wings of 
silicate emission unless the sources are OH/IR stars.) The color temperature Td is calculated by 
constructing a look-up table with the excess ratio calculated as per Equation [5] for a wide range of 
temperatures, and selecting the value of temperature from this table that reproduces the observed 
ratio of excesses. Once the color temperature is calculated. Equation [3] can be solved for the optical 
depth. The color temperature and optical depth values calculated using these simplified equations 
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are primarily useful in investigating the trends of color temperature and optical depth with excess. 
The equations break down for sources with high optical depth (i.e., for extreme AGB stars). The 
temperature derived from a ratio of two broad band fluxes is a simplification that will in turn affect 
the optical depth calculation. We will model the CS shells of our AGB candidates considering 
details such as the wavelength dependence of the emissivity and the choice of dust species in a 
future paper to obtain more precise estimates for the temperatures and optical depths. 

Since we are only interested in the overall trend of the optical depth with excess flux, we 
fix Rin at ten times the stellar radius i?*, which is calculated by scaling the radii of the model 
photospheres to the luminosity of each sta r. According t o models of dust condensation, most of the 



dust has condensed by ~ lOi?* (see, e.g.. iHofner 1120071 ). In practice, the dust condensation radii 



and Rin/R* values will vary, but this complication is ignored in this work as we are only interested 
in a comparative study of the MIR color temperatures and optical depths. While the temperature 
of CS dust varies with radius from ~1000 to ~100 K, the (single) color temperature calculated 
here will be dominated by the distance corresponding to the hottest and most optically thick dust 
emission in the CS shell, which is typically close to the inner radius of the dust shell. Our values 
for describe the optical depth of the high density inner region of the CS shell. This estimate is 
insensitive to large but cool optical depths. This optical depth will be the dominant component of 
the line of sight optical depth. Table [3] shows the color temperatures and optical depths for the 
sources in Tables [1] and [2l Data for the entire set of sources with valid 8 yum excesses is available 
in the electronic version of Table [3l 



4. Results 

In this section, we investigate the relation of the MIR excess to other observed parameters for the 
sample including luminosity, color temperature of the CS dust, and opacity of the CS dust shell. 
We focus on the statistical nature of the sample and global trends of the results. This approach 
provides the range of physical parameters for the sample and is a first and necessary step to obtain 
appropriate ranges for model parameters for any detailed modeling that may follow for the sample. 



4.1. Excess-Luminosity Relations 

The calculated IR excesses show a general increasing trend with luminosity at all IRAC wavelengths 
for O- and C-rich AGB sources. This trend is most obvious at 8 fim (Figures [7] and [8]) . While 
the correlation is similar at the shorter IRAC wavelengths, the number of sources with significant 
reliable excesses is considerably smaller. This is expected for the warmer, moderately obscured 
non-extreme sources. The 8 and 24 /xm excesses of the extreme AGB candidates also correlate well 
with their luminosities, as is evident from Figure [H There is a huge spread in 24 /xm excesses for 
the 0-rich and C-rich candidates overall, but an increasing trend is apparent for the 0-rich AGB 
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candidates at luminosities above ~ 3 x 10'* Lq. The large spread at low luminosities may be caused 
by significant variations in the MLRs during the early stages of AGB stars. The plots of 0-rich 
sources show the bright and faint population sources in light and dark grey respectively. We also 
display the median SED for each population as an inset. For the C-rich and extreme sources, we 
show the median SED for three equally-populated bins in luminosity. Relations of the form 

logX,(mJy) = A + 51ogL(L0) (6) 

are fit to the 8 and 24 /im excesses, resulting in the following best-fit parameters: 

A = -6.7 ± 0.5, B = l.7± 0.1 (O - rich, 8 ^m) 
A = -9.4 ± 3.2, B = 2.l± 0.7 (O - rich, 24 ^m, L > 3 x 10* L©) 
A = -6.6 ± 6.2, B = 1.7 ± 1.5 (C - rich, 8 ^m) 
A = -7.0 ±15, B = 1.6 ± 3.5 (C - rich, 24 ^m) 
A = -5.0 ± 0.3, S = 1.5 ± 0.1 (Extreme, 8 /um) 
A = -8.3 ± 0.7, B = 2A± 0.2 (Extreme, 24 /um) 

The poor fit to the 24 ^m excess for the C-rich candidates reflects the huge spread of 24 fim 
excesses for these sources. 



4.2. Color Temperature and Opacity 

Figures [10] and [11] show the variation of the color temperature and 24 /xm optical depth with 24 
^m excess for all three types of AGB candidates. As expected from their redder [8] — [24] color, 
the members of the faint 0-rich population tend to have cooler temperatures and higher optical 
depths than the bright population, and they also show a more pronounced variation with excess. 
There is a huge spread in the C-rich AGB color temperatures. The maximum color temperature is 
higher than in the 0-rich case, suggesting that the carbonaceous dust grains become hotter than 
the silicates. AGB winds composed of oxy gen-rich compounds are less efficient at absorbing visible 



photons than those that are carbon-rich (jWallerstein fc Knapplll998l ). Carbonaceous grains are 



efficient absorbers of optical photons and are highly emissive at IR wavelengths in comparison to 
silicates, thus they are more likely to reach higher temperatures. We fit the color temperatures for 
the extreme AGB candidates with a power-law relation of the form 

logTdiK) = A + BlogXUmSy) 

and obtain the following best-fit values: 

A = 2.9± 0.4, B = -0.19 ± 0.02 (Extreme) 

while a similar relation, 

log r24^m = A + BlogX^ (mJy) (7) 
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when fit to the 24 fj,m optical depths, gives 

A = -0.49 ± 0.12, B = 1.6± 0.1 (O - rich, faint) 
A = -1.2 ± 0.3, B = 0.9± 0.2 (O - rich, bright) 
A = -1.7 ± 0.3, B = 1.2± 0.2 (C - rich) 
A = -1.8 ± 0.1, S = 0.79 ± 0.09 (Extreme) 

Figures [12] and [13] show the variation of the 24 /im optical depth with luminosity and color temper- 
ature respectively. There is no trend in the optical depth of the 0~rich candidates with luminosity. 
For the C-rich sources, there is no correlation, but there appears to be an absence of sources 
with low optical depth at higher luminosities. For the extreme AGB stars, there appears to be a 
positive correlation of higher optical depth with higher luminosities. The optical depth decreases 
monotonically with the color temperature for all three types of AGB candidates. CS dust shells 
have temperature gradients from ~ 1000 K in the interior to ^ 100 K in the outer regions. With 
increasing optical depth, therefore, we are probing the progressively cooler, outer regions of the 
AGB star candidates. 



Discussion 



The 8 and 24 /xm band excess fluxes show an increase with increasing luminosity of the central 
star. For each type of AGB candidate, the slope of the excess-luminosity relation is similar across 
the four IRAC and MIPS24 bands, though in the IRAC 8 /im and MIPS24 bands, excesses are 
"reliable" only beyond ~0.1 mJy (see Figure [6]). The increase of exc ess with luminosity is consistent 
with the MLR-luminosity rela tion found bvlvan Loon et al.l (| 19991 ). The mass loss-rate is roughly 
proportionate to tL (see, e.g., Ilvezic fc Elitzui]ll99d ^. so we expect that the MLR increases with 
increasing luminosity as long as the optical depth does not decrease faster than ~L^^. This suggests 
that the excess is a good reflection of the MLR. More quantitative comparisons of MLR for this 
whole sample with luminosity is beyond the scope of this study, but this will be addressed in a 
future paper. 

The uncertainties in photometry alone cannot account for the considerable spread in the calcu- 
lated excesses. Stars on the AGB suffer from variability and episodic mass loss, and our observations 
of these stars may be capturing their fluxes at different epochs in their variability cycle, or during 
episodes of increased or diminished mass loss. At a given luni i nosity , these effects introduce a 
variation in the MLR and hence also affect the excess. IVijh et al.l (120091 ) show a median variability 
index of about 3 for the 0-rich and C-rich AGB sources and about 5 for the extreme sources. This 
would correspond to a change in the excess fluxes of about the same factor. The color temperature 
depends on the ratio of the excesses, and the effect of variability is probably weaker. The effects of 
grain chemistry, metallicity of the environment, and the C/0 ratio of these sources have not been 
accounted for in this study, as this requires spectroscopic determination. At a given luminosity, 
there is also a dependence on progenitor mass - a degeneracy between a less evolved, more massive 
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star and a more evolved, less massive star arises. Such degeneracy effects increase the spread at 
the lower luminosity end (only the most massive stars evolve to high luminosities). 



Observations for Galactic AGB stars (jGuandalini et al.ll2006l : iBusso et al.ll2007l ) show that C- 
rich AGB stars are in general more obscured than their 0-rich counterparts, but this is partly due 
to the higher opacity of carbonaceous dust. Note that extreme 0-rich AGB stars {i.e., OH/IR 
stars) are even more obscured than these C-rich sources. In the case of the extreme sources, there 
is a high optical depth regardless of chemistry of the CSE. As an AGB star evolves, more mass in 
the form of dust is deposited into the surrounding shell, increasing its opacity and causing further 
reddening of stellar light and stronger emission at longer wavelengths, corresponding to cooler 
color temperatures. The color temperature and optical depth variations for extreme AGB stars in 
Figures [10] and [TT] support these claims. The existence of a redder, fainter population of 0-rich 
sources first seen in the IRAC CMD is also apparent in plots of the excess, color temperature 
and optical depth. The faint population is in general cooler and more obscured than the bright 
population. The SAGE survey is able to differentiate between these two types of sources for the first 
time. The bright population corresponds to the young, most massive AGB stars that prevent the 
dredge-up of carbon by undergoing hot-bottom burning (part of seq uence G inlNikolaev &: Weinberg 
([20001)). This population also emerges in the model isochrones of [Marigo et al.[ ([2008[ ) - they find 
that the thermally-pulsing (TP) phase in their isochrone of age log(i/yr) = 8.2 is well-developed 
and populated by 0-rich stars undergoing HBB. The faint population consists of AGB stars in 
which HBB does not occur, res ulting in a smooth transition fro m 0-rich to C-rich ch e mistr y. 
The log(t/yr) = 9.1 isochrone of [Marigo et al.[ ([200a) (sequence J in [Nikolaev &: Weinberg! (j200d )) 
represents the typical evolutionary phase corresponding to this population. 

Our fj-clipping method of recovering sources with reliable excesses can be used t o compare our 



resul ts with the expected lower-limit to measurement of MLRs by the SAGE survey ([Meixner et al 
2006[ ). The mid- infrared CMDs in Figures [H[ [T5l and [T6[ show the sources with reliable 8 fiui 



excesses as black circles. We detect 0-rich and C-rich sources with reliable excesses up to within 
^0.1 mag of t he tip of the RGB ([8.0]~11.9). which is below the detection limit mentioned in 
Meixner et al.l ([2006 ) ([8.0]=11.0) for measuring significant mass loss. The extreme AGB stars in 
our sample are all well above this detection limit, showing that the SAGE survey detects all the 
extreme mass- losing AGB sources in the LMC. 

We est i mate the current LMC AGB st ar mass-loss bu d get u sing a method analogous to 
Blum et al.l ([20071 ). The MLRs obtained by [van Loon et al.[ ([l999l ) for spectroscopically identi- 
fied AGB stars are plotted against their SAGE 8 fim excess fiuxes in Figure [T7[ The figure also 
shows power-law fits of the form 



logM(10-^ Meyr" 



A + BlogX^{mJy) 



to all three types of sources with the following best-fit parameters: 

A = -1.5, B = 1.0 (0-rich) 
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A = -1.7, B = 1.1 (C - rich) 
A = -2.8, B = 1.7 (Extreme) 

These relations are then used to derive MLRs for ah the sources in our Usts with reUable excesses. 
Figure [18] shows the total dust MLR as a function of LMC AGB luminosity. The extreme AGB 
stars, despite their comparatively low numbers, are the most significant contributors to the AGB 
mass-loss budget in the LMC. We find the dust injection rates from 0-rich, C-rich and extreme 
AGB stars to be 0.14, 0.24, and 2.36 xlO~^ Mq yr~^ respectively, which puts the estimate for the 
total dust injection rate from AGB stars at 2.74 x 10~^ Mq yr~^. While stars with luminosities 
between M;,o; = —7.1 and M;,o; —7.8 may be hi ghly embedded AG B stars, sources brighter than 
Mboi > —7.8 are definitely supergiants (see, e.g.. ISloan et al.ll2009l ). The luminosity distribution 
(Figure H]) shows that there are very few sources in this range, but it is clear from Figure [18] that 
the contribution of the red supergiants to the dust injection rate is comparable to that of the AGB 
stars. This introduces a significant uncertainty in the value for the AGB dust injection rate derived 
in this section. 

Calculation of the total (gas+dust) injection rate requires knowledge of the gas:dust ratio, 'i/'j 
in the circumstellar shells of AGB stars. The lower metallicity of the LMC will result in a higher 
silicate gasidust ratio but the carbon gas:dust ratio may be similar to Galactic values (IHabing 
Keeping this in mind, we choose gasidust ratios of 200 an d 500 for C-ri c h and 0-rich 



stars respectivelj0. The C-rich value is similar to that obtained by Skinner et al. ( 1999 ) for the 
Galactic carbon star, IRC+10 216. The gas MLR from the 0-rich and C-rich AGB candidates 
is then 7xl0~^ Mq yr~^ and 4.8x10"'^ Mq yr~^ respectively. The chemical identification of the 
extreme AGB sources is not possible from our data, but assuming ip=200 and ■0=500 for these 
stars will provide lower and upper limits. The gas MLR for the extreme AGB population is in 
the range 4.7-11. 8x 10"^ yj.-!. Thus, the total AGB gas MLR is (5.9-13)xl0"3 yr-i. 
We will improve on our estimate with ou r follow-up radiate tra nsfer modeling of the dust shells 
around these stars using the 2DUST code (|Ueta &: Meixneiil2003l ). The results of the present study 
will help constrain the range of optical depths and grain temperatures. The optical depth is an 
important input parameter for the code. 

Taking into account the effects of variability on the excesses as discussed previously, we estimate 
that the MLR of an individual source can vary on average by a factor of ~3 for the 0-rich and 
C-rich sources and ~15 for the extreme AGB candidates. Due to the large number of 0-rich and 
C-rich sources, the errors of the flux variations on the cumulative MLR estimates would probably 
cancel out. For the extreme AGB candidates, the larger effects of the variability and the smaller 
numbers may result in a significant change in the cumulative MLR. The magnitude of such an effect 
could be calibrated with an IR variability study of the LMC AGB stars to determine their periods 
and hence phase-correct the fluxes and excesses, but no such study exists for our complete sample. 



"The work of 



van Loon et al. ( 2008f ) suggests that ip scales in proportion to the metallicity for both 0-rich and 



C-rich AGB stars, but this effect can not be accounted for in our simple estimate. 



- 15 - 



Our C-rich AGB star gas MLR o f 4.8x10 M0 yr is cor nparable to the value, 6 x10 Mg 





yr~^, found by iMatsuura et al.l (j2009l ). However, we note that iMatsuura et al.l (|2009l ) base their 
measurement on MLR versus color relations and include the extreme C-rich AGB stars, identified 
by IRS spectroscopy, as well as the color classified C-rich AGB stars. Our excess versus MLR 
relation method is two orders of ma gnitude more sen s itive to lower MLR sources (dust MLRs of 
~ 10"" Mq yi-^) compared to the IMatsuura et al.l ^2Qod ) method (dust MLRs of 10"^ M© 
yr"^). It is therefore fortuitous that our estimates are so similar and our exclusion of the highest 
MLR extreme C-rich AGB stars is balanced by the inclusion of more low MLR C-rich AGB stars. 
Furthermore, the total mass injection rate of C-rich mass loss is probably higher than both these 
estimates, but definitive numbers will require a clear identification of what fraction of the extreme 
AGB stars are C-rich. 

Our estimate for the total L MC AGB gas MLR is higher than those obtained for two lo- 



cal group dwarf irreg ulars, WLM (jjackson et al 
(|.Tackson et allboOTbl (0.2-1.0) x 10"^ M© yr^^). 



2007a 



(0.7-2.4) x 10-3 M0 yr"^) and IC 1613 
As their values for the MLRs of individual stars 



are in good agreement with the estimates of individual stars by Ivan Loon et al.l (jl999l ) for the LMC, 
the disparity between our results and the values for WLM and IC 1613 could arise due to various 
reasons. The total MLR return due to AGB stars depend s on the total number of A GB stars in a 
galaxy. The total stellar mass of the LMC 3 x 10^ M^, van der Marel et al. 2002) is two or ders 



of magnitude greater than those estimated by iJackson et al.l (j2007al ) and iJackson et al.l (|2007bl ) for 
WLM and IC 1 613 (1.1 x 10^ Mp, and 1.7 x 10^ Mq respectively). Moreover, the age of the LMC 
bar is 4-6 Gyr (jSmecker-Hane et al.ll2002l ). which is optimal for an enhanced AGB population at 



presen t times, compared to the slightly older stellar population of WLM (9 Gyr, IJackson et al 



2007al ). Based on the differences in stellar masses and ages, we expect that the LMC has consider- 
ably more AGB stars than either WLM or IC 1316, leading to a higher total MLR from AGB stars. 
However, the differences in the MLR determination may also increase o ur MLR compared to th eirs. 
Our IR excess method is able to detect sources with small MLR, while Ivan Loon et al.l (|l999l ) and 
Jackson et al.l (j2007al ) derive their estimates from the highest mass-losing sources alone. In this 
sen se, ours is a more comp lete estimate of the total MLR. On the other hand, our interpolation of 
the Ivan Loon et al.l (|l999l ) MLRs may not be appropriate for sources with very low excess. Thus, 
while the excess method is able to identify sources with low MLRs, the derived MLRs may be 
overestimates. 

The current star formation rate (SFR) of the LMC is estimated to be about 0.1 Mq yr~^ 
(jWhitnev et al.l l2008l ) . which is an order of magnitude higher than our calculated AGB MLR. 
Thus, the current star formation rate is not sustainable, assuming the LMC is a closed box system, 
and that AGB stars are the only means of replenishing the ISM. However, the most massive stars 
{e.g., red supergiants, supernova explosions) may also contribute to the mass-loss return, and these 
need to be accounted for to complete the census of the mass budget return to the LMC. In addition, 
interactions with the Small Magellanic Cloud may cause infall or loss of mass that would need to 
be considered in the mass budget of the ISM. 



6. Summary and conclusions 



We classify evolved stars in the SAGE Epoch 1 Archive based on their 2MASS and IRAC colors, 
and use photospheric models to estimate the infrared excess emission due to circumstellar dust. 
We obtain about 8200 0-rich, 5800 C-rich, and 1400 extreme AGB sources with reliable 8 jum 
excesses (SNR > 3). The corresponding numbers in the 24 //m band are about 4700, 4900, and 
1300 respectively. The excesses increase with an increase in luminosity in all four IRAC bands as 
well as the MIPS24 band. We use the 8 and 24 /xm excess fluxes to derive dust color temperatures 
and optical depths, and we observe that higher excess fluxes correspond to cooler temperatures 
and optically thicker dust shells. These quantities for our list of AGB candidates with valid 8 
fim excesses are available in the form of electronic tables. We also estimate the present day AGB 
mass-loss budget in the LMC by comparing modeled mass-loss rates with our excesses estimates 
to find that the extreme AGB stars are the most significant contributors to mass loss in the LMC. 
Our data also suggests that the rate of dust injection from red supergiants to the LMC ISM is 
comparable to the total dust injection rate from the AGB stars, which we calculate to be about 
(5.9-13) X 10-^ M0 yr-^ 

The authors wish to thank Rita Gautschy-Loidl and Peter Hauschildt for providing us with 
their photospheric models. We would also like to thank the anonymous referee for their helpful 
comments. The SAGE Project is supported by NASA/Spitzer grant 1275598 and NASA NAG5- 
12595. 
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Table 1. Source List of Color-Selected AGB Star Candidates ^ 



Identifier'' 






Typeb 


magU'^ 


(5magU 


ma 


gB 


5ma 


gB 


magV 


<5ma 


gV 


magI 


(5magl 


magj 


(5magJ 


ma 


gH 


<5ma 


gH 


magK 


(5magK 


SSTISAGEIA 


J054938 


72- 


-683458 


2 


O-rich 


99.99 


99.99 


19 


46 





04 


17.55 





05 


13.76 





07 


12 


01 


0.03 


11 


11 





03 


10 


78 


0.02 


SSTISAGEIA 


J055530 


35- 


-684647 


5 


O-rich 


20.46 


0.19 


18 


11 





03 


16.18 





03 


13.95 





04 


12 


59 


0.02 


11 


72 





02 


11 


48 


0.02 


SSTISAGEIA 


J055420 


11- 


-680449 


5 


O-rich 


99.99 


99.99 


18 


87 





04 


16.59 





05 


13.53 





03 


11 


91 


0.02 


11 


02 





02 


10 


71 


0.03 


SSTISAGEIA 


J055729 


20- 


-684444 


2 


O-rich 


17.97 


0.06 


17 


47 





04 


16.32 





08 


16.38 





30 


11 


93 


0.02 


11 


03 





02 


10 


75 


0.02 


SSTISAGEIA 


J055321 


17- 


-683114 


7 


O-rich 


21.04 


0.24 


18 


49 





09 


16.51 





05 


13.87 





03 


12 


52 


0.02 


11 


60 





02 


11 


37 


0.02 


SSTISAGEIA 


J054522 


57- 


-684244 


5 


C-rich 


99.99 


99.99 


20 


68 





06 


16.45 





03 


13.82 





04 


12 


49 


0.02 


11 


17 





02 


10 


38 


0.02 


SSTISAGEIA 


J055650 


80- 


-675030 


5 


C-rich 


99.99 


99.99 


20 


42 





07 


16.66 





11 


13.56 





04 


12 


12 


0.02 


11 


01 





03 


10 


35 


0.02 


SSTISAGEIA 


J055835 


31- 


-682009 


7 


C-rich 


99.99 


99.99 


20 


13 





05 


16.28 





03 


13.55 





04 


12 


30 


0.02 


11 


02 





03 


10 


20 


0.02 


SSTISAGEIA 


J055311 


71- 


-684720 


9 


C-rich 


99.99 


99.99 


21 


91 





17 


17.93 





04 


14.56 





04 


12 


03 


0.02 


10 


91 





03 


10 


14 


0.02 


SSTISAGEIA 


J055036 


67- 


-682852 


3 


C-rich 


21.61 


0.44 


19 


31 





05 


16.68 





03 


14.42 





04 


12 


87 


0.03 


11 


92 





03 


11 


52 


0.02 


SSTISAGEIA 


J052742 


48- 


-695251 


5 


Extreme 


19.08 


0.09 


18 


58 





09 


18.72 





11 


14.86 





07 


14 


13 


0.07 


12 


41 





06 


11 


08 


0.04 


SSTISAGEIA 


J052714 


19- 


-695524 


3 


Extreme 


20.04 


0.12 


19 


65 





07 


18.36 





05 


14.45 





05 


13 


74 


0.03 


12 


13 





03 


10 


85 


0.02 


SSTISAGEIA 


J053239 


06- 


-700157 


5 


Extreme 


99.99 


99.99 


20 


29 





20 


17.30 





10 


14.81 





04 


13 


16 


0.03 


11 


73 





02 


10 


73 


0.02 


SSTISAGEIA 


J052950 


52- 


-700000 


1 


Extreme 


18.80 


0.09 


18 


28 





09 


17.40 





14 


16.08 





05 


13 


96 


0.04 


12 


30 





04 


10 


80 


0.03 


SSTISAGEIA 


J053441 


38- 


-692630 


7 


Extreme 


99.99 


99.99 


22 


78 





45 


20.26 





11 


14.59 





04 


12 


36 


0.02 


10 


92 





02 


9 


87 


0.03 



^The electronic version of this tabic contains all the AGB candidates considered in this study. 
^SAGE Epoch 1 Archive (SAGEIA) designation, including position coordinates of IRAC source. 

''Sources are classified as O-rich or C-rich based on their J-Kg colors, or as "Extreme" based on their 2MASS and IRAC colors. For more details, see 32] 
'^Magnitudes and errors in the UBVI, JHKs , IRAC and MIPS 24 ^m bands. A value of 99.99 in any band represents either a saturation or a non-detection. 



Table 1. (Continued) 



mag36 


i5mag36 


mag45 


(5mag45 


mag58 


i5mag58 


mag80 


5mag80 


mag24 


i5mag24 


10.53 


0.03 


10.72 


0.03 


10.50 


0.04 


10.44 


0.05 


9.86 


0.07 


11.34 


0.03 


11.40 


0.03 


11.28 


0.04 


11.19 


0.04 


10.59 


0.14 


10.47 


0.04 


10.64 


0.02 


10.46 


0.05 


10.40 


0.04 


9.97 


0.07 


10.47 


0.03 


10.60 


0.02 


10.39 


0.03 


10.30 


0.04 


10.01 


0.08 


11.07 


0.07 


11.02 


0.04 


10.90 


0.03 


10.76 


0.03 


10.10 


0.08 


9.67 


0.03 


9.70 


0.03 


9.68 


0.04 


9.30 


0.04 


8.89 


0.04 


9.76 


0.04 


9.86 


0.03 


9.81 


0.04 


9.30 


0.04 


9.26 


0.03 


9.34 


0.05 


9.03 


0.02 


8.80 


0.03 


8.68 


0.03 


8.61 


0.03 


9.51 


0.05 


9.43 


0.04 


9.27 


0.04 


8.99 


0.03 


8.66 


0.03 


11.06 


0.04 


11.08 


0.04 


10.89 


0.04 


10.68 


0.05 


10.21 


0.09 


9.02 


0.05 


8.65 


0.04 


8.39 


0.03 


8.00 


0.03 


7.49 


0.02 


9.62 


0.04 


9.12 


0.03 


8.71 


0.03 


8.24 


0.03 


7.61 


0.03 


9.38 


0.04 


9.33 


0.03 


9.19 


0.04 


8.83 


0.03 


8.46 


0.03 


8.94 


0.04 


8.23 


0.04 


7.55 


0.03 


6.85 


0.02 


6.03 


0.02 


9.06 


0.05 


8.33 


0.05 


7.81 


0.02 


7.29 


0.03 


6.74 


0.02 



Tabic 2. Mid-infrared Excess FlTixes 



Identifier Type X36^ SX36^ X45 5X45 X58 5X58 X80 5X80 X24 5X24 

(mjy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) 



SSTISAGEIA 


J054938 


72-683458 


2 


O— rich 


4 


546 





622 


1 


918 





315 


2 


338 





287 


1 


335 





210 





457 





051 


SSTISAGEIA 


J055530 


35-684647 


5 


O— rich 





843 





287 





727 





171 





726 





135 





450 





091 





213 





053 


SSTISAGEIA 


J055420 


11-680449 


5 


O-rich 


4 


254 





784 


1 


932 





279 


2 


170 





363 


1 


222 





187 





349 





050 


SSTISAGEIA 


J055729 


20-684444 


2 


O-rich 


4 


351 





531 


2 


413 





292 


2 


740 





230 


1 


693 





190 





326 





050 


SSTISAGEIA 


J055321 


17-683114 


7 


O-rich 


2 


387 





654 


2 


392 





268 


1 


909 





163 


1 


313 





100 





428 





049 


SSTISAGEIA 


J054522 


57-684244 


5 


C-rich 


18 


980 


1 


265 


12 


830 





720 


8 


010 





533 


7 


247 





488 


1 


283 





078 


SSTISAGEIA 


J055650 


80-675030 


5 


C-rich 


13 


170 


1 


375 


7 


974 





645 


5 


065 





509 


6 


534 





504 





606 





048 


SSTISAGEIA 


J055835 


31-682009 


7 


C-rich 


29 


860 


2 


606 


31 


540 





989 


26 


200 





957 


16 


000 





617 


1 


761 





076 


SSTISAGEIA 


J055311 


71-684720 


9 


C-rich 


19 


900 


1 


948 


16 


800 


1 


066 


13 


170 





819 


9 


984 





475 


1 


570 





077 


SSTISAGEIA 


J055036 


67-682852 


3 


C-rich 


1 


080 





494 


1 


209 





257 


1 


308 





207 





953 





160 





238 





051 


SSTISAGEIA 


J052742 


48-695251 


5 


Extreme 


69 


590 


2 


958 


62 


220 


2 


082 


50 


750 


1 


422 


40 


600 


1 


094 


7 


258 





160 


SSTISAGEIA 


J052714 


19-695524 


3 


Extreme 


39 


680 


1 


490 


40 


240 


1 


117 


37 


570 


1 


141 


32 


440 





895 


6 


457 





154 


SSTISAGEIA 


J053239 


06-700157 


5 


Extreme 


49 


920 


1 


958 


33 


430 





812 


24 


320 





834 


18 


790 





459 


2 


950 





068 


SSTISAGEIA 


J052950 


52-700000 


1 


Extreme 


74 


320 


3 


098 


91 


780 


3 


533 


110 


200 


2 


910 


117 


100 


2 


549 


27 


820 





396 


SSTISAGEIA 


J053441 


38-692630 


7 


Extreme 


66 


460 


3 


197 


83 


610 


3 


607 


86 


330 


1 


828 


77 


850 


2 


089 


14 


490 





305 



^The electronic version of this table contains all the AGB candidates with valid 8 /nm excesses. 
''The MIPS 24 /nm and IRAC 8.0, 5.8, 4.5, and 3.6 /xm band excess fiuxes and errors in mjy. 
''The errors in the excess fluxes have been calculated by propagating the photometric errors: 
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Table 3. Color temperatures and Optical Depths 



Identifier 




Type 


(K) 


ST 
(K) 


r(8 tJ,m)^ 


<5t(8 /um) 


t(24 /im) 


<5t(24 /im) 


SSTISAGEIA 


J054938, 


.72-683458.2 


O-rich 


393 


28 


0.813 


0.017 


0.938 


0.003 


SSTISAGEIA 


J055530, 


.35-684647.5 


O-rich 


352 


37 


0.806 


0.028 


0.935 


0.005 


SSTISAGEIA 


J055420, 


.11-680449.5 


O— rich 


422 


36 


0.885 


0.012 


0.962 


0.002 


SSTISAGEIA 


J055729, 


.20-684444.2 


O-rich 


502 


48 


0.916 


0.008 


0.972 


0.002 


SSTISAGEIA 


J055321, 


.17-683114.7 


O-rich 


401 


21 


0.743 


0.016 


0.914 


0.003 


SSTISAGEIA 


J054522, 


.57-684244.5 


C-rich 


525 


26 


0.922 


0.004 


0.974 


0.001 


SSTISAGEIA 


J055650, 


.80-675030.5 


C-rich 


814 


84 


0.982 


0.001 


0.994 


0.000 


SSTISAGEIA 


J055835, 


.31-682009.7 


C-rich 


705 


32 


0.953 


0.002 


0.985 


0.000 


SSTISAGEIA 


J055311, 


.71-684720.9 


C-rich 


559 


22 


0.931 


0.002 


0.977 


0.000 


SSTISAGEIA 


J055036 


.67-682852.3 


C-rich 


446 


53 


0.949 


0.007 


0.983 


0.001 


SSTISAGEIA 


J052742, 


.48-695251.5 


Extreme 


521 


10 


0.830 


0.003 


0.944 


0.001 


SSTISAGEIA 


J052714 


.19-695524.3 


Extreme 


495 


9 


0.793 


0.004 


0.931 


0.001 


SSTISAGEIA 


J053239, 


.06-700157.5 


Extreme 


559 


11 


0.908 


0.002 


0.970 


0.000 


SSTISAGEIA 


J052950, 


.52-700000.1 


Extreme 


456 


5 


0.601 


0.005 


0.867 


0.001 


SSTISAGEIA 


J053441, 


.38-692630.7 


Extreme 


511 


9 


0.808 


0.003 


0.936 


0.001 



^The electronic version of this table contains all the AGB candidates with valid 8 //m excesses. 

'''Color temperatures and related uncertainties, derived from the 8 /nm and 24 /im excesses. The uncertainties are 
calculated by propagating the photometric errors. 

''Optical depths and related uncertainties at 24 /im and 8 /tm derived from the color temperature. The uncer- 
tainties are calculated by propagating the photometric errors. 



1.0 1.5 2.0 2.5 3.0 

J-Ks 



Fig. 1. — A Hess diagram showing the locus of AGB stars in the Kg versus J-Kg CMD. Stars 
brighter than the tip of the RGB (Kg =12, dashed hne) are select ed as AGB candid ates, and these 



are classified as 0-rich and C-rich based on the CMD cuts in ICioni et al.l (|2006l ). The O-rich 



sources are bounded by the lines Ki and K2 (shown as solid line in figure) in their paper, while 
sources redward of K2 are classified as C-rich. The extreme AGB candidates (not shown in this 
figure) are classified based on their 2MASS and IRAC colors. 
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• O-rlch AGB 

• C-rich AGB 

■ Extreme AGB 
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1 2 
[3.6]-[8.0] 



2. — The three types of AGB candidates (0-rich: blue, C-rich: red, Extreme: green) on an 
vs [3.6]-[8.0] CMD. The tip of the RGB is at [3.6]Ri[8.0]=11.9. 
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Fig. 3. — The 0-rich (blue), C-rich (red) and extreme AGB (green) stars on the [8.0] -[2 4] versus [24] 



CMP . The sources below the solid line are the fainter, redder population mentioned in lBlum et al 



(|2006l ) (The finger "F" in their Figure 6). About 80% of the 0-rich stars in our sample belong to 
this population. 
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Fig. 4. — The luminosity function for our sample of AGB candidates (blue: 0-rich, red: C-rich, 
green: Extreme, black: total). The vertical dashed line is the classical AGB luminosity limit. 
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Fig. 5. — The panels on the top show the normahzed SEDs for the best-fit 0-rich (left) and C-rich 
(right) models in light gray. The response curves for the optical UVBI, 2MASS JHKg , IRAC and 
MIPS24 filters are also shown. Folding the transmission profiles of each filter into the model SED 
results in a flux in every fllter (red diamonds). The bottom panels show the median flux in each 
band for the hundred bluest (in V-Kg color) 0-rich (left) and C-rich (right) sources (black squares) 
plotted over the flux of the corresponding best-fit model (red diamonds) scaled to the median flux 
in the H band. The H band magnitudes of these bluest sources are in the range 11.99 to 11.36 
(0-rich) and 11.95 to 10.32 (C-rich). The 0-rich model is a 1 model with r=4000 K and 
logff = 0, while the C-rich model has r=3200 K, a C/0 ratio of 1.3, and log^ = -0.43. 
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Fig. 6. — Histogram of the 8 /xm excesses for the faintest 0-rich AGB candidates (and thus with 
the lowest S/N ratio excess emission). The solid line is for all sources. The dotted line shows 
sources with S/N>3. The dashed line shows sources with excesses considered unreliable and hence 
were excluded from our study. 
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10^ 10* 10^ 10^ 

Lumlnoalty (L^ 

Fig. 7. — The 8 nm (top) and 24 (bottom) excesses plotted against luminosity for the O- 
rich sources. Sources from the bright and faint populations (light gray and dark gray circles 
respectively) are shown. The solid lines show, respectively, power-law fits to all the sources with 
reliable 8 fim excesses (top) and sources with reliable 24 //m excesses brighter than 3 x 10^ Lq 
(bottom). Representative error bars are also shown. Inset: the median SEDs of the bright and 
faint populations. 
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io* 10^ 
Lumlnoalty (L^ 

Fig. 8. — The 8 //m (top) and 24 nm (bottom) excesses plotted against luminosity for the C-rich 
sources. The solid lines are power-law fits to the excess-luminosity relations. (The fit is very poor 
in the 24 fim case) Representative error bars are also shown. Inset: the median SEDs for three 
equally-populated luminosity bins (solid line: least luminous one-third of the sample, thin dashed 
line: intermediate luminosity, thick dashed line: most luminous). 
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Fig. 9. — The 8 /um (top) and 24 //m (bottom) excesses plotted against luminosity for the extreme 
sources. The solid lines are power-law fits to the excess-luminosity relations. Representative error 
bars are also shown. Inset: the median SEDs for three equally-populated luminosity bins (solid 
line: least luminous third of sample, thin dashed line: intermediate luminosity, thick dashed line: 
most luminous). 
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24 fim Excess (mjy) 

Fig. 10. — The color temperature derived from the ratio of the 8 /xm and 24 /im excesses plotted 
against the 24 /xm excess for all three types of AGB candidates (blue: faint 0-rich, cyan: bright 
0-rich, red: C-rich, green: Extreme). The cross is a representative error bar. 
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Fig. 11. — The 24 fim optical depth derived from the 8 and 24 excesses for all three types 
of AGB candidates. The color coding is the same as for Figure [TOl Power-law fits to the optical 
depth-excess relation are shown (dotted: faint 0-rich, dashed: bright 0-rich, dot-dashed: C-rich, 
solid: Extreme). The cross is a representative error bar. 



- 35 - 




Fig. 12. — The 24 /um optical depth as a function of luminosity for all three types of AGB candidates. 
The color coding is the same as for Figure [TOl 
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Fig. 13. — The 24 //m optical depth plotted against color temperature for all three AGB types. 
The color coding is the same as for Figure [101 



t I I I U I I I u 

1 2 3 1 2 3 
[3,6]-[8,0] [8,0]-[24] 

Fig. 14.— [8.0] vs [3.6]-[8.0] and [8.0] vs [8.0]-[24] CMDs for the 0-rich AGB stars in our list 
with 8 and 24 /xm detections (gray circles). The sources with "reliable" excesses are superimposed 
(black circles). The tip of the RGB is at [8.0]=11.9 (dashed line) 
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Fig. 15. — Same as Figure [14] for our C-rich AGB candidates 
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Fig. 16. — Same as Figure [14] for our extreme AGB candidates 
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8 fim Excess (mjy) 



Fig. 17. — The mass- loss rates from Ivan Loon et al.l (|l999l ) for 0-rich (triangles), C-rich (filled 
squares) and extreme (circles) AGB sources plotted against their 8 fim excesses. The lines show 
power-law fits for the dust mass-loss rate-excess relation for all three types of AGB candidates 
(solid line: 0-rich, dashed line: C-rich, dot-dashed line: Extreme). 
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10^ lO'* 10* 10® 

Luminosity (L©) 



Fig. 18. — The cumulative dust mass-loss rate as a function of luminosity for the 0-rich (thin solid 
line), C-rich (dashed line) and Extreme (dot-dashed line) AGB stars in our lists. The thick solid 
line is the total AGB dust mass-loss rate as a function of luminosity. The extreme AGB stars are 
the major contributors to the mass-loss rate, at 2.36x10^^ Mq yr~^. The total dust mass-loss rate 
is 2.74x10"^ Mq yr~^. The vertical dotted line is the classical AGB luminosity limit. As is evident 
from the figure, there are very few, if any, carbon-rich AGB candidates above this limit. However, 
deeply embedded extreme AGB stars and 0-rich stars undergoing hot-bottom burning can exceed 
this luminosity limit. 



